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We measure the branching fraction for the charmless semi-inclusive process B r/'Xs, where 
the 77' meson has a momentum in the range 2.0 to 2.7 GeV/c in the T(4S) ccnter-of-mass frame 
and Xs represents a system comprising a kaon and zero to four pions. We find B{B — > rj'Xs) = 
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(3.9 ± 0.8(stat) ± 0.5(syst) ± 0.8(model)) x 10"*. We also obtain the Xs mass distribution and find 
that it tends to favor models predicting high masses. 



The production of high momentum rj' mesons in B me- 
son decays is expected to be dominated by the B t]'Xs 
process, where Xs is a strange hadronic system, gen- 
erated by the b sg* transition as depicted in Fig. 
n^a-c). Figure n^d) shows the color-suppressed modes 
_^ rj'D^*^^, which are significant sources of back- 
ground and which have been measured for the first time 
recently 0. Contributions from b —> u transitions and 
other sources of rj' are expected to be negligible 0- 




u,d - 



FIG. 1: Lowest order diagrams for (a,b,c) B r)' Xs and (d) 
the color-suppressed background B" ri'D^*^'\ 

The large B — > tj' Xs branching fraction measured by 
the CLEO collaboration |2j, prompted intense theoret- 
ical activity, which focused the special character of the 
rj' meson as receiving much of its mass from the QCD 
anomaly . 

A later measurement by CLEO confirmed the large 
ij' production, measuring B{B rj'Xnc) = (4.6 ± 
l.l(stat) ± 0.4(syst) ± 0.5(bkg)) x 10"^ 8], where X^c 
denotes a charmless recoiling hadronic system. 
The rate for B — > rj'Xs and especially the fully 
background-subtracted distribution of the mass of Xs 
can provide important clues to the dynamics of weak de- 
cays and to the structure of the isosinglet pseudoscalar 
mesons. 

We present results for the branching fraction B{B ^ 
rj'Xs) and the mass spectrum of Xs- The signal is ana- 
lyzed for rj' momentumbetween 2.0 and 2.7 GeV/c in the 
CM to suppress background coming from b ^ c rj' 
cascades such as i? ^ DsX with Dg v'^i B DX 
with D rj'X, B A^X with Ac rj' X. Our analysis 
is based on data collected with the BABMi detector Q at 
the PEP-II asymmetric e+e~ collider located at the Stan- 
ford Linear Accelerator Center. An integrated luminosity 
of 81.4 fb~^, corresponding to 88.4 million BB pairs, was 
recorded at the T(4S') resonance (on-resonance) and 9.6 



fb^^ were recorded 40 MeV below this resonance (off- 
resonance), for continuum background studies. 

Two tracking devices are used for the detection of 
charged particles: a silicon vertex tracker consisting of 
five layers of double-sided silicon microstrip detectors, 
and a 40-layer central drift chamber, both operating in 
the 1.5 T magnetic field of a superconducting solenoid. 
Photons and electrons are detected by a CsI(Tl) electro- 
magnetic calorimeter. Charged-particle identification is 
provided by the average energy loss (dE/dx) in the track- 
ing devices, and by an internally refiecting ring-imaging 
Cherenkov detector covering the central region. 

We select BB events by requiring at least four charged 
tracks and a value of the ratio of the second to zeroth 
Fox- Wolfram moment ^3 less than 0.5. 

We form a B candidate by combining an ?/ — > ?77r+7r~, 
where the rj decays into 77, with a K'^ or a that is 
reconstructed in the tt+tt" channel, and up to four pi- 
ous, of which at most one is a 7r°, leading to 16 possible 
channels 



B" 

B+ ^Tj'K°TT+i+Tr°) B" 
B+ rj'K°7r+T:+7r-{+n") B" - 



r;'if07r+7r-(+7r") 



The mass of the 77 —^ 77, K'^ — > tt+tt^ and tt^ —^ 77 
candidates are required to lie within 3ct (ct = 16, 3 and 
6 MeV/c^ respectively) of their known values and are 
then kinematically constrained to their nominal masses. 

To identify the s quark in the Xg system, we re- 
quire a -fCg or a track consistent with a charged kaon. 
The charged-kaon selection has been optimized to reduce 
background from B — s- rfn, rj' p, and rj' ai decays. For the 
ifg, we require the angle a between the momentum of the 

candidate and its fiight direction to be less than 0.05 
radians, as it peaks at zero for true if^ particles. 

We require candidates for B rj'Xg to be consistent 
with a B decay, based on the beam-energy-substituted 
mass, mES = \/ {s/2 + pa.pB)^/E^ - p| and the energy 
difference, Ai? — Eg — ^fs/l, where E and p denote the 
energy and momentum of the particles, the subscripts 
and B refer to the initial T{AS) and the B candidate, re- 
spectively, the asterisk denotes the T(4S') rest frame, and 

is the e^e~ center-of-mass energy. In addition, the 
cosine of the angle between the thrust axis of the B candi- 
date and that of the rest of the event in the center-of-mass 
frame (cosfl^) is used to remove continuum background, 
which is peaked near | cos 6*^ | = 1 , while signal events are 
uniformly distributed. We require mES > 5.265 GeV/c^, 
|Ai?| < 0.1 GeV, and | cos6'5'| < 0.8. For each event, we 
select the candidate with the smallest , with defined 
by 
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m(T|ira) (GeV/c^) m(rii7t) (GeV/c^) 

FIG. 2: Fits to the r/Trvr invariant mass for on-resonance (top) 
and off-resonance (bottom) data samples, for the modes (a,b) 
/s:± and (c,d) ifg. 

- (mEs - Ms)V^'(™es) + {^Ef/aH^E) , 

where is the i?-nieson mass and where where the res- 
olutions (7(mEs) = 3 MeV/c^ and o-(A£') = 25 MeV are 
obtained from Monte Carlo simulation. The remaining 
continuum background is subtracted with the use of off- 
resonance data. 

The background contribution from color-suppressed 
modes rj'D^*^^ is estimated from a Monte Carlo 

simulation which uses our measurement of its branching 
fraction, ^ V^^*^") ^ (1.7±0.4(stat)±0.2(syst)) x 

10-4 1]. 

To determine efficiencies, we model the signal using a 
combination of the two-body mode B — > rj' K and, for Xg 
masses above the Kir thresholdja non-resonant derived 
from the theoretical predictions maiEl, which are based 
on the anomalous ?7'-gluon-gluon coupling and which fa- 
vor high-mass Xs systems. The fraction of the two-body 
mode is constrained in the simulation model to be be- 
tween 10% and 15% [H E^. When not forming a K 
meson, the X^ fragments into sq and sqg {q = u,d). We 
find that the overall efficiency is (6.0 ± 0.2)% for the 
modes and (4.7 ± 0.1)% for the modes, including the 
branching fraction B{K^ tt+tt"). 

The branching fraction of _B ^ rj'Xg is computed 
through a fit to the number of if signal events, with 
rj' momentum between 2.0 and 2.7 GeV/c, both for on- 
resonance and off-resonance data. To parameterize the 
background, we use a Gaussian function for the signal 
and a second order polynomial. For the fit of the off- 
resonance data sample, we constrain the mass and width 
of the rj' to the values obtained with on-resonance data. 
Figure 121 shows the fits of the rjinr invariant mass distri- 
butions for the and K'^ modes. The fitted yields are 
reported in Table |l| 
The semi-inclusive branching fraction is computed by 



TABLE I: Resuhs of the fits for and ifl? modes. Yields 
for on-resonance data (Von), off-resonance data (Yoff), ex- 
pectation from color-suppressed background (ybs) and on- 
resonance data after background subtraction (Y) are given. 
A luminosity scale factor, / — 8.48 , is applied to the off- 
resonance yield. 





modes 


/('s modes 


Yon 


577.0 ± 34.0 


367.0 ± 34.0 


VOFF 


18.9 ±8.5 


21.7 ±8.4 


Yes 


63.6 ± 11.4 


26.9 ±4.5 


Y 


353.1 ±80.5 


156.1 ±79.1 



TABLE 11: Contribution of different sources to the systematic 
error for modes with a K"^ or K^. 



Source 


syst (%) K'i 


syst (%) 


Tracking 


3.4 


3.3 


ri, tt" detection 


7.0 


8.2 


K/Kl ID 


2.5 


4.3 


3.4 


3.4 


Nbb 


1.1 


1.1 


MC sample size 


3.0 


3.0 


jy'£)(*)o subtraction 


3.0 


2.9 


Total 


12.1 


13.5 


Model 


20 


20 



performing a weighted average of the results obtained 
for the and modes. The detection efficiencies 
are corrected to account for the r/ and branching frac- 
tions to the channel we observe. For the modes, we 
convert the result so it corresponds to and K . The 
final state Xs includes both K^- and if°-tagged decays. 
Assuming that their branching fractions are equal, we 
obtain B{B rj'X,) = (3.9 ± 0.8(stat) ± 0.5(syst) ± 
0.8(model)) x 10"''. We obtain the systematic error by 
combining the sources listed in Table Hll 

The largest uncertainty arises from our model of the 
Xg system. To estimate that uncertainty, we use an al- 
ternative model which consists of a combination of res- 
onant modes: rj'K, rj'K* (892), r]'Ki{1270), ?7'i4:i(1400), 
r7'ii:*(1410), 77'^2*(1430), r]'K*{1780), and 7fK*{204:5). 
The variability of the efficiency and our knowledge of 
the resonant sector lead us to assign a 20% system- 
atic uncertainty. Other systematic uncertainties include 
track reconstruction efficiency, reconstruction efficiencies 
of 7r° — > 77, r] 77, and K'^ — > 7r"'"7r~ candidates, 
charged-kaon identification efficiency, secondary branch- 
ing fractions, number of BB events (Ngg), the size of our 
Monte-Carlo sample, and subtraction of the background 
from B° r7'D(*)o. 

To explore the Xs mass distribution, we select B can- 
didates for which the mass of the 77' is within three stan- 
dard deviations of the known value, and subtract the 
continuum contribution by using on-resonance data in 
the sideband 5.200 < tues < 5.265 GeV/c^. The contin- 
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(a) 



s 

o 20 - 




0.5 I 1.5 
m(Xj) (GeV/c^) 



(b) 


If 







m(X.) (GeV/c-) 



FIG. 3: Continuum-subtracted K rni invariant-mass distribu- 
tions for (a) all B modes and (b) B'^ modes, including combi- 
natorial background. Solid and dashed histograms represent 
expected backgrounds from r]'D° and rj'D*^, 

respectively. 



■ K'(892) m K',(1430) 

i K,( 1270)0 K',(1780) 

T K,(1400)D K'.(2045) 

* K'(1410) 



2.5 



0.5 I 1.5 2 
m(X ) (GeV/c-) 



2.5 



FIG. 4: Variation of efficiency with m{Xs). The filled circles 
indicate the efficiency for non-resonant Xs simulation. The 
other symbols denote the values for the resonances. 



uum background scaling factor {A), from the sideband 
to signal regions, is computed from off-resonance data 
to be 0.591 ± 0.118. The resulting mass distributions are 
shown in Fig. |31for all B modes and separately for the 
modes. The peak at m{Xs) — 500MeV/c^ corresponds 
to the two body mode B rj'K. 

To obtain the full Xs spectrum, we fit the rj' mass 
distribution in bins of Xs mass. The efficiency, aver- 
aged over the charged and neutral kaons, as a function 
of ■m{Xs), is shown in Fig. ^ The correction for the 
feed-across between bins is included in the efficiencies. 

According to simulations, the Xs system is correctly 
reconstructed for 85% (60%) of the candidates in the re- 
gion miXs) < 1.5 GeV/c2 (m(X^) > 1.5 GeV/c^). For 
correctly reconstructed events, the experimental resolu- 
tion varies from 5 to 15 MeV/c^ for low and high masses, 
respectively. In the case of misreconstructed events, the 
resolution ranges from 100 to 150 MeV/c^. Table Hill 
shows the fitted yields for the raw signal, the sideband 
region, the expected color-suppressed background, and 
the yield after full background subtraction, as a function 



TABLE III: Fitted yields for on-resonance data and color- 
suppressed background for different m{Xs) ranges in GeV/c^. 
The sideband yields (Ysb) must be corrected by the sideband 
to signal region scaling factor (see text) before subtraction. 



m{Xs) ranj 


;e Yon 


Ysb 


Yes 


Y 


[0.4,0.6] 


200 ± 15 


46.1 ±8.8 




172.8 ±15.9 


[0.6, 1.2] 


120 ± 14 


100 ± 13 




60.9 ± 16.0 


[1.2,1.5] 


114 ± 15 


112 ± 14 


1.1 ±0.3 


46.7 ± 17.1 


[1.5,1.8] 


150 ± 18 


163 ± 17 


7.7 ± 1.6 


46.0 ± 20.7 


[1.8,2.0] 


140 ± 17 


93± 15 


47.4 ± 9.6 


37.6 ± 21.4 


[2.0,2.3] 


149 ± 20 


142 ± 18 


26.2 ±4.5 


38.9 ±23.1 


[2.3,2.5] 


80 ± 14 


70± 14 


4.9 ±0.9 


33.7 ± 16.3 




(b) 




m(X J (GeV/c-) 



m(X.) (GeV/c") 



FIG. 5: Branching fractions as a function of m{Xs). Both 

(a) and (b) show the same data, though the efficiency used 
in (a) is derived from the non-resonant model, while that in 

(b) the efficiency comes from the model with a combination 
of resonances. The errors include bin-to-bin systematics; an 
additional systematic error of ~ 8% (not shown) is common 
to all points, (a) The open histogram represents the expec- 
tation from non-resonant m{Xs) simulation, (b) The open 
histogram represents the expectation from a mixture of reso- 
nant modes with equal proportions. The hatched histogram 
results if some heavy resonances are enhanced. The equal 
mixture provides a good approximation to what is predicted 
in JJ. 



of m{Xs). 

The branching fraction as a function of 'm(Xs), ob- 
tained from the fully background-subtracted yield (Table 
IIII|I . is shown in Fig. 

We compare data and simulation by forming a 
difference. The probability for the nonresonant Xg 
model (Fig. E^a)) to fit the data is 61% while it is close 
to ~ 10~^ for the equal mixture of resonances (Fig. [Sfb)). 
We find improved agreement with the resonant model if 
the weights of and are increased by a factor of 
1.5, leading to a probability of 2%. 

As a consistency check of the method, we measure the 
two-body decay modes {Xs = K'^jK^), and find 171. 0± 
14.0 and 27.1 ± 5.6 events in on-resonance data for rj'K^ 
and jy'i^s respectively, and no rj' signal events for both 
channels in off-resonance data, leading to the branching 
fractions B{B^ -> r]'K'^) = (6.9 ± 0.6(stat)) x lO''^ and 
B{B° f]'K°) = (5.6± 1.2(stat)) x 10"^ These values 
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are fully compatible with what has been measured by re- 
cent exclusive analyses 0, . 

In summary, we have measured the branching fraction, 
B{B ~* r]'Xs) = (3.9±0.8(stat)±0.5(syst)±0.8(model))x 
10-^ for 2.0 < p*{r]') < 2.7 GeV/c. We have also de- 
rived the m{Xs) spectrum and found that the data tends 
to confirm models predicting a peak at high masses and 
seems to disfavor predictions based only on the diagram 
of Fig. ma,b) for which m{Xs) peaks near 1.4-1.5 GeV/c^ 

M 

Among the various theoretical conjectures to explain 
this production, an Tj'gg coupling due to the QCD 
anomaly has been widely suggested as a likely explana- 
tion. However, the r]'gg form factor initially proposed Q 
is disfavored by recent studies of the inclusive production 
r{lS) r]'X [ll[l3. A recently updated approach ^ 
exploiting the same tj' gluon anomaly could in princi- 
ple account for the observed branching fraction and the 
m{Xs) spectrum. 
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